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Abstract — An attempt has been made to study the unsteady 
MHD two-dimensional free convective elastico-viscous fluid past 
an inclined permeable plate with Dufour effects in presence of 
chemical reaction. The elastico-viscous fluid flow is 
characterized by Walters liquid (Model B'). The perturbation 
scheme is used to solve the governing equations of the fluid 
motion. The approximate solutions have been derived for the 
velocity profile, temperature field, concentration field, skin 
friction, rate of heat transfer and rate of mass transfer. It is 
noticed that the flow field is considerably affected by the 
elastico-viscous parameter. The velocity profile and the 
approximate skin friction coefficient have been presented 
graphically to observe the elastico-viscous effects for various 
values of the flow parameters involved in the solution. 


Index Terms — Chemical Reaction, Dufour effects, MHD, 
Skin friction 


I. INTRODUCTION 

The analysis of free convection flow adjacent to inclined 
surface bounded by an extensive body of fluid is of 
considerable importance in micrometeorology and industrial 
applications. Investigation of chemical reaction, heat and 
mass transfer in MHD elastico-viscous fluid flow for an 
electrically conducting fluid past a heated surface has 
attracted the interest of many researchers in view of its 
applications in many engineering problems such as plasma 
studies, petroleum industries, MHD power generator, cooling 
of nuclear reactors, MHD propulsion thermo fluid dynamics, 
boundary layer control in aerodynamics, etc. Process 
involving heat and mass transfer in porous media are often 
encountered in the chemical industry and in reservoir 
engineering in connection with thermal recovery process. 

A comprehensive review of the studies of heat transfer 
mechanism through porous media has been made by Neild 
and Bejan [1]. The study of natural convection flow of 
inclined plate has been presented by the authors Ganesan and 
Palani [2] and Sparrow and Husar [3]. Free convection 
boundary layer flow over a horizontal and slightly inclined 
surface has been studied by Pera and Gebhart [4]. Karkus [5] 
has applied the perturbation technique to study the natural 
convection flow adjacent to inclined isothermal and 
finite-length surfaces. Hossain et al. [6] have investigated the 
free convection flow from an isothermal plate inclined at a 
small angle to the horizontal. Raptis et al [7] have investigated 
the viscous flow over a non-linearly stretching sheet in 
presence of chemical reaction and magnetic field. Rajesh et al 
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[8] have studied the chemical reaction and radiation effects on 
MHD flow past an infinite vertical plate with variable 
temperature. Diffusion thermo and radiation effects on MHD 
free convective heat and mass transfer flow past an infinite 
plate in the presence of chemical reaction of first order has 
been studied by Raveendra Babu et al [9]. The heat and Mass 
Transfer in MHD Boundary Layer Flow past an inclined plate 
with viscous dissipation in porous medium has been analyzed 
by Singh P K [10]. Considering the model of elastico-viscous 
fluid, many scientists have solved problems of engineering 
interests viz. Saxena and Dubey [11], Choudhury and 
Mahanta [12],Choudhury and Dey [13], Choudhury and Dhar 
[14]. 

The aim of the present work is to study the diffusion thermo 
(Dufour effect) and radiative effects on MHD unsteady free 
convective flow of elastico-viscous fluid past an inclined 
permeable plate characterized by Walters liquid (Model B') . 
The effects of the elastic-viscous fluid on velocity component 
and skin friction coefficient have been shown graphically with 
the combination of other flow parameters involved in the 
solution. 

The constitutive equation for Walters liquid (Model B') is 
u* = -pg ik + 2ij 0 e LI1 - 2k v e r<k (1) 


where f 1 * is the stress tensor, p is isotropic pressure, S\ k is 
the metric tensor of a fixed co-ordinate systemr ! , ir is the 
velocity vector, the contravariant form of 9 "^ is given by 


9 ~ St 


? rr, ik 

' ,TJl 


( 2 ) 


It is the convected derivative of the deformation rate 

tensor 9 [hl defined by 

2e !K = (3) 


Here % is the limiting viscosity at the small rate of shear 
which is given by 

JJ-u = Jjj N (t 1 dr and = J 0 tN CV) dr, (4) 

iV (V) being the relaxation spectrum as introduced by Walters 

[15]-[16]. This idealized model is a valid approximation of 

Walters liquid (Model B') taking very short memories into 

account so that terms involving 

_f D T n N(:r)dT, n> 2 (5) 


have been neglected. 

The mixture of Polymethyl Methacrylate and pyridine at 250 
C containing 30.5 gm of polymer per litre and having density 
0.98 gm/ml. fits very nearly to this model. For this mixture, 
the relaxation spectrum as given by Walters is 

NCX) = a%S<X) +■ —^% (0 < A < p] 

= 0 T < |3 

where cr = 0.13,% =7.9 poises (gm/cm. sec), |3 = 0.18 sec 

and is the Dirac’s delta function so that k0=0.60557 

gm/cm. 
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II. Mathematical Formulation 

The region of unsteady two-dimensional free convective 
hydromagnetic flow of an elastico-viscous fluid past an 
infinite inclined permeable plate has been analyzed in 
presence of chemical reaction. A uniform magnetic field B- is 
applied normal to the plate. Here the mechanism of heat and 
mass transfer in connection with Dufour effect are studied. 
Since the plate is of infinite length, all physical variables are 
functions of y and t only. Let x r -axis is taken along the 
inclined plate and y '-axis is taken normal to the plate. Let u r 
and v r be the velocity components in x* and y r direction 
respectively. Let a be the small angle made by the plate with 
the vertical so that sina=0 as shown in Figure 1. 



where G r is the Grashof number for heat transfer, G m is the 
Grashof number for mass transfer, Bn is the magnetic field, 
8 is the volumetric co-efficient of expansion for heat transfer, 
8 is the volumetric co-efficient of expansion for mass transfer, 
g is the acceleration due to gravity, T and T = are the 
temperature of the fluid inside the thermal boundary layer and 
the fluid temperature in the free stream respectively, c and c = 
are the corresponding concentrations, IV, is the co-efficient 
of mass diffusivity, A T r is chemical reaction parameter, K is 
the permeability parameter. 

Introducing the non-dimensional parameters in the equations 
(7), ( 8 ), (9), we get 

flu du d 2 u / dh 0*u\ / \\ 

s _L ^ =+ ^ _ki W- L ^) He " B+ ^“ B -r + ir) u (12) 

30 30 13 Z 0 3 Z 0 

3t ~ L lty = '^6y* +D ' 1 3^ (13 ^ 

30 30 1 3 : 0 

3t L 3y S c dy 2 

where k ± is the visco-elastic parameter. 


(14) 


Figure 1: Geometry of the problem 


The boundary conditions of the problem in the dimensionless 
form are 

u = 0.. 0 = 1, 0 = 1 at y = 0 

u = 0.. 0=0, 0 = 0 as y — do (15) 


Under the above assumptions the flow is governed by the 
following set of equations: 


Equation of continuity: 
dv 

— = 0 
dy 

v = constant = —i? E , 
Equation of motion: 


du _ du d 2 u jfcfl 
di V dy V dy 2 p 

—- u 

P 




d d u _ d d u 

didy 2 11 dy 3 


( 6 ) 

F ,g8(;T - fj) c o s 

(7) 


Energy equation: 

dT dT 

T7fv- = 


dT dT k d 2 T IV d 2 C 


dt ' ' dy pCpdy 2 + C s C p dy 2 

Concentration equation: 
dc dc d 2 c 

With relative boundary conditions: 
u = IT, v = v(t), c = cL aty = 0 
u = 0, T = T.,., c = cL as y -> do 
The non-dimensional parameters_are^_ 


( 8 ) 


(9) 


( 10 ) 
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a/?v(T w -TJ gpvic^-cj) _uCp 
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n v D m K r ( c^-cj uB$v 

S c = 7 =-—--M = 




KU 2 vK l k 0 U 2 v, 


p U 2 ' 


K = ~^r,K r = ~rrr,k = ^r,L = 77 

U ■* 


u 2 


pv £ 


( 11 ) 


III. METHOD OF SOLUTION 

We assume the solutions of the non-linear partial differential 
equations (12) to (14) of the following form 
u = u 0 (y) F ftiUifyX 0 = & 0 iy} F (y) 

0 = 0o (y) +*i0i (16) 

Using (16) in the equations (12), (13) and (14) and equating 
the coefficient of like powers of k ± and neglecting the highest 
powers of we get 

ug^ -G r %cosa - G m <p D cosa (17) 
u[ + Lu ‘- (m F — juj. = -G r 9 L cosa - V 0 iCdsg - (18) 


KJ 

F = — 

8 [ + P r L8l= -P r D u <p[ 

0 jj "1“ LS C 0 p — S € K r 0 □ = 0 
0i F LS c <p[ - 5 c A r r 0i = 0 
The modified boundary conditions are 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 


ii 0 = 1 , — 1 , 0 c = l f u 1 = 0 , = O, 0 L = 0 at y = 0 

%, = 0,00 = O,0 E = Q,-^ = 0,^ = 13,0! = 0 as y-t k (23) 
Solving equations (21) and (22 ) and using boundary 
conditions (23), we get 

0* (24) 

0i = 0 (25) 

Solving equations (19) and (20) and using boundary 
conditions ( 23), we get 

0 O =A 1 e- a & -A 2 e-^' (26) 

9 t = 0 ^ (27) 

Solving equations (17) and (18) and using boundary 
conditions (23), we get 

u 2 = F^p"^' (28) 

iq = Aje - ^ 2 y — A^e~^ sy F A L ^e~ asy (29) 

The constant are obtained but not presented here for the sake 
of brevity. 
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IV. RESULTS AND DISCUSSION 


The velocity profile for the fluid flow is given by 

u - 4 { 4 e_ffsJ + +■ ) (30) 


The skin friction coefficient at the plate in non-dimensional 
form is given by 

[ du id fdu\ dim 

(31) 

The non-dimensional form of the rate of heat transfer 
coefficient in terms of Nusselt number is given by 
1 fd&\ 1 fd% d& L \ 

t32) 

The non-dimensional form of the rate of mass transfer 
coefficient in terms of Sherwood number is given by 


1 fdtp\ 1 fd<p 0 dtp t ' 

S e V dy Kl dy 


- _ (i!t) 
~S c \dy) ir=0 ~ 


(33) 


y=o 


The purpose of this study is to bring out the effects of 
elastic-viscous parameter on the unsteady MHD free 
convective flow past an inclined permeable plate with Dufour 
effects in presence of chemical reaction. The elastico-viscous 
effects is exhibited through the non-dimensional parameter 
.The corresponding results for Newtonian fluid are 
obtained by setting = 0 . 

To understand the physics of the problem the fluid velocity u 
is depicted against y in Fig. 1 to Fig. 6 for different values of 
Prandtl number P r , Magnetic number M, Grashof number for 
heat transfer G r , Grashof number for mass number G m , Dufour 
number D u , Chemical reaction K r , Schmidt number Sc and 
elastico-viscous parameter ki. The numerical calculations are 
carried out for Permeability parameter K=2 and L=2. Fig. 1 
and Fig. 2 reveals that an increase in the Grashof number does 
not bring any significant change in the flow pattern but with 
the rising effect of visco-elasticity the fluid velocity enhances 
in comparison with Newtonian case. It is also observed that 
the fluid velocity diminishes due to the rise of magnetic 
parameter (Fig. 1 and Fig. 3) and Dufour number (Fig. 1 and 
Fig. 4) in both Newtonian and non-Newtonian cases. The 
variation of chemical reaction (Fig. 1 and Fig. 5) and Schmidt 
number (Fig. 1 and Fig. 6) shows that the fluid velocity 
enhances in both Newtonian and non-Newtonian cases. 

Fig. 7 to Fig. 12 depict the variation of the shearing stress for 
different values of the flow parameters involved in the 
solution. Fig. 7 reveals as Prandtl number increases viscous 
drag declines for both Newtonian and non-Newtonian 
parameter. It is observed from Fig. 8 that the increasing values 
of Magnetic parameter amplify magnitude of shearing stress. 
Fig. 9 illustrates that the rising values of Dufour number 
increase the shearing stress for both Newtonian and 
non-Newtonian fluid flow phenomenon. The variation of 
shearing stress against chemical reaction K r (Fig. 10) shows 
that the growth of visco-elastic parameter subdues the 
shearing stress. Fig. 11 shows that the magnitude of shearing 
stress increases with amplified values of Grashof number for 
cooling plate (G r >0) but it decreases with the increasing 
values of elastico-viscous parameter in comparison with 
Newtonian fluid.. The reverse pattern of flow is observed for 
the heated plate (G r <0) (Fig. 12). It is further noticed that the 
rate of heat and mass transfer coefficients are not significantly 
affected by the elastico-viscous parameter. 


V. CONCLUSIONS 

The unsteady two-dimensional free convective 
elastico-viscous fluid past an inclined permeable plate with 
Dufour effects in presence of chemical reaction has been 
investigated for different values of flow parameters. In the 
analysis, the following conclusions are made: 

■ The fluid velocity diminishes due to the rise of magnetic 
parameter and Dufour number in both Newtonian and 
non-Newtonian cases. 

■ The speed of both Newtonian and non-Newtonian fluids 
through the inclined permeable plate enhances with the 
increasing values of chemical reaction and Schmidt number. 

■ The growth of Grashof number does not bring any 
significant change in the flow pattern of the elastic-viscous 
fluid. 

■ Viscous drag declines for both Newtonian and 
non-Newtonian fluid with the increase of Prandtl number and 
chemical reaction parameter. 

■ The increasing values of Magnetic parameter amplify 
magnitude of shearing stress. 

■ The rising values of Dufour number increase the shearing 
stress for both Newtonian and non-Newtonian fluid. 

■ The shearing stress enhances with amplified values of 
Grashof number for cooled and heated plate. 



Fig: 1. Variation of u against y for P r =5, G r =6, G m =7, M = 3, 
K = L= 2, D u =0.03, K r = 0.6, S c =3 



Fig: 2. Variation of u against y for P r =5, G r =10, G m =7, M = 3, 
K = L= 2, D u =0.03, K r = 0.6, S c =3. 
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Fig: 3. Variation of u against y for P r =5, G r =6, G m =7, M = 5, 
K = L= 2, D u =0.03, K r = 0.6, S c =3 


Fig: 7. Variation of x against P r for G r = 3, G m =8, M = 3, K 
=L=2, D u =0.03, K r =0.6, S c =3. 



Fig: 4. Variation of u against y for P r =5, G r =6, G m =7, M = 3, 
K = L= 2, D u =0.05, K r = 0.6, S c =3. 



Variation of x against M for P r =5, G r = 3, G m =8, K=L=2, 
D u =0.03, K r =0.6, S c =3. 



Fig: 5. Variation of u against y for P r =5, G r =6, G m =7, M = 3, 
K = L= 2, D u =0.03, K r = 0.9, S c =3. 



Fig: 9. Variation of x against D u for P r =5, G r = 3, G m =8, M=3, 
K=L=2, K r =0.6, S c =3. 
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Fig: 6. Variation of u against y tor P r =5, G r =6, G m =9, M = 3, 
K = L= 2, D u =0.03, K r = 0.6, S c =4. 


Fig: 10. Variation of x against K r for P r =5, G r = 3, G m =8, M = 
3, K =L=2, D u =0.03, S c =3. 
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Fig: 11. Variation of x against G r (> 0) for P r =5, G m =8, M = 3, 
K =L=2, D u =0.03, K r =0.6, S c =3. 
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Fig: 12. Variation of x against G r (< 0)for P r =5, G m =8, M = 3, 
K =L=2, D u =0.03, K r =0.6, S c =3. 
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